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Abstract 
CrGeTezrecently emerges as a new two-dimensional (2D) ferromagnetic 
semiconductorthat is promising for spintronic device applications. Unlike 
CrSiTe;whose magnetism can be understood using the 2D-Ising model, 
CrGeTe3exhibits a smaller van der Waals gap and larger cleavage energy,which could 
lead to a transition of magnetic mechanismfrom 2Dto 3D.To confirm this speculation, 
we investigate the critical behaviorCrGeTe; around the second-order 
paramagnetic-ferromagnetic phase transition. We obtain the critical exponents 
estimated by several commonexperimental techniques including the modified Arrott 
plot, Kouvel-Fisher method and critical isotherm analysis, which showthat the 
magnetism of CrGeTesfollows the tricritical mean-field model with the critical 
exponents p, y, and ô of 0.240 + 0.006, 1.000 + 0.005, and5.070 + 0.006, respectively, 
at the Curie temperature of 67.9 K. We therefore suggest that the magnetic phase 
transition from 2D to 3D for CrGeTesshould locatenear a tricritical point. Our 
experiment providesa direct demonstration of the applicability of the tricritical 


mean-field model to a 2D ferromagnetic semiconductor. 


Since the successful exfoliation of single layer graphene, 2D materials have been 
attracting significant interest due to their highly tunable physical properties and 
immense potential in scalable device applications [1-5]. However, pristine graphene 
exhibits no bandgap and its inherent inversion symmetry suppresses the spin-orbit 
coupling (SOC)[6-11]. The weak SOC and zero bandgapeliminate graphene as a 
potential candidate forbeing applied to spintronic devices,which require one to search 
for alternative 2D materials that extend beyond graphene to other layered 
materialswith van der Waals gaps [7-11]. For example, in single-layer MoSz, the large 
SOC leads to a unique spin-valley coupling which may be useful for spintronic 
applications [12-16]. Whereas spintronic devicesusing 2D materials arestill in their 
infancy[17-20], which is due to the lack of long-range ferromagnetic order that is 
crucial for macroscopic magnetic effects[21,22].The emergence of ferromagnetism in 
2D materials in combinationwith their rich electrical and optical properties could open 
up ample opportunities for 2D magnetic, magneto-electric, and magneto-optic 
applications[18,19,23]. 

Recently,Chromium TelluridesCrXTe; (X = Si, Ge, and Sn) with the 
centrosymmetric have arisen significant attention because they belong to a rare 
category of ferromagnetic semiconductors possessing a 2D layered structure [7,23-35]. 
Extensive theoretical and experiment efforts have been extended towards 
understanding the properties of these 2D magnets. On the theoretical side, recent 
studies on CrXTe3 have been focusing on their electronic structure and magnetic 


properties, particularly predictions of the single-layer properties[24,28-33]. On the 


experimental side, CrSiTe; and CrGeTe; have been prepared and 
characterized[7,23,25-27,34,35]. Comparing with CrSiTe3, showing characteristics of 
a 2D-Ising behavior[7,34,35], the smaller van der Waals gap and the larger in-plane 
nearest-neighbor Cr-Cr distance in CrGeTe; enhance the Curie temperature from 32 K 
for the CrSiTe; to 61 K for the CrGeTe3[7,25,28,31]. In addition, theoretical 
investigations have suggested that the single-layerCrGeTe; presents characteristics of 
2D-Ising behavior similar to CrSiTe3[31,33]. By contrast,a scanning magneto-optic 
Kerr microscopy experiment, single-layerCrGeTe; represents a close-to-ideal 2D 
Heisenberg ferromagnetic system using the rigorous renormalized spin wave theory 
analysis and calculations[23]. It is known that, withthe increase of the X atom radius, 
CrXTe3 presents the smaller van der Waals gap and the larger cleavage energy 
[7,25,28,31]. We supposethat CrXTe3 system may undergo a three dimensional (3D) 
magnetic phase transition from 2D with the increase of the X atom radius. Therefore, 
a method to rapidly characterize thecritical behaviorof single-crystalline CrGeTesis 
crucial.For this purpose, we present a detailed investigation of the critical phenomena 
of CrGeTe3using the initial isothermal M(H) curves around the Curie temperature Tc. 
We find that the critical exponents of CrGeTe3 satisfy the universality class of the 
tricritical mean-field theory. This indicates that the magnetic phase transition of 
CrGeTe; should be close to a tricritical point from 2D to 3D. 

Samples of single-crystalline CrGeTe; were prepared by the self-flux 
technique[26].The XRDdata indicated that the powders are singlephase with the 


rhombohedral structure (see Supporting Information). We measured the heat capacity 


using the Quantum Design physical properties measurement system (PPMS-9T) and 
characterized the magnetic properties by the magnetic property measurement system 
(MPMS-XL5).Density functional theory (DFT) calculations were performed using 
Vienna Ab-initio Simulation Package[36]. We used the local density 
approximation[37,38] to treat the electron-electron exchange-correlation interactions. 
The electron-ion interactions are described by the potentials based on the projector 
augmented wave method[39,40]. 

Figure 1 (a) and (b) show the inverse of temperature-dependent susceptibility 
1/y(T) of CrGeTe3 under field cooled cooling with applied magnetic field H = 100 
Oe, parallel to the ab plane and c axis, respectively. We observe a 
paramagnetic—ferromagnetic (PM-FM) transition that occurs at a critical temperature 
T? of 67.3 K,as determined by the derivative of the susceptibility. This temperature is 
consistent with the values of 61 K or 70 K reported previously[25-27]. For a FM 
system, the1/y (T) above T."* can be described by the Curie-Weiss law resulting from 
the mean-fieldtheory[41]. The red curves showing the Curie-Weiss law are obeyed 
only at high temperature. A close observation of Fig.1 (a) and (b) reveals that the 
1/ x(T) curves deviate from straight lines at around 150 K, which is much higher than 
T."* , indicating strong short-range FM spin interactions in CrGeTezabove T."* . The 
effective magnetic moment is determined to be H, ~ 4.22 44, (parallel the ab plane) 
and 4.35 44, (parallel the c axis), which are close to the theoretical value expected for 
Cr” of 3.87 Hp, [25]. The insets of Fig.1 (a) and (b)show the isothermal 


magnetizationM(H) at 5 K exhibiting a typical FM behavior with the saturation field 


Hsof about 5 kOe (parallel to the ab plane) and 2.5 kOe (parallel to the c axis). In 
addition, the M(H)curves show almost no coercive force for CrGeTe3. 
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Figure 1.(Color online): (a) and (b) Inverse of temperature-dependent susceptibility 1/7 (T) of 


CrGeTe3 under field cooled cooling with an applied magnetic field H of 100 Oe, parallel to the ab 
plane and c axis, respectively. The red solid lines are the fitted results according to the Curie-Weiss law. 


The insets show the isothermal magnetization curvesM(H) at 5 K.(c) Specific heat C, as a function of T 
for CrGeTe3 and the fitted on rove (T) using Eq. (1) and (2); Temperature-dependent magnetic (d) 


specific heat Cmag(T) and (e) entropy Smag(T). The blue dashed line refers to Smag(T— œ) calculated with 


the magnetic moment S of Cr?” being 3/2. 

Figure 1(c) shows the variation of the zero-field specific heat (SH)C,(T) with 
temperature. The sharp anomaly in C,(T) at 64.8 K corresponds to the Curie 
temperature T .Since CrGeTe3 is a seminconductor[25],the electronic contribution to 
the heat capacity is not considered. The Cmagcan be calculated by the following 


equations[41]: 


Cms (T) = C,(T)- NC?™ (T), (1) 
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where R is the molar gas constant,@pis the Debye temperature, andN = 5 is the 
number of atoms per formula unit. The sum of Debye functions accounts for the 
lattice contribution to specific heat. We can extract the magnetic contribution Cmag(T) 
from the measured specific heat of CrGeTe3. The fitted C,(T) for CrGeTe3 by Eq. 
(1)and (2) over the temperature range from about 7 to 250 K is shown by the red 
curve in Fig.1(c) using the Debye temperatureOp= 476.5 K. We observe a sharp peak 
at T" of 64.8 K and there is strong dynamic short-range FM spin interactionsabove 
TS" (shown in the Fig.1(d)). The magnetic entropy Smag(T)is calculated by 


C (T 
Smag (T) = N Cos (P) gy .Fig.1 (e) shows the temperature dependence of Smag(T). The 


T 

entropy of CrGeTe3;per mole with completely disordered spins S is 
S mag (T > ©) = 6RIn(2S+1) .Using S = 3/2 for Cr**, we obtain Smag(T — œ) of 69.2 
J/(mol-K).As shown in Fig.1(e), one observes that the 64.7 J/(mol-K) of Smag(T) at 150 
K is close to the value ofSmag(T — œ), which suggests that CrGeTe3 presents the 
obvious PM behavior above 150 K. However, we observe the Smag is 36.3 J/(mol‘K) at 
T2", which is only 52.5% of the value for Smag(T — œ) in Fig.1(e). These results 
indicate there is strong dynamic short-range FM spin interactions between T$" and 
150 K, which is consistent with the above conclusion from fitting the1/y(T ) data to 
the Curie-Weiss law. 

As mentioned above, with the X atom radiusincreases,the CrXTe3compounds 


present the smaller van der Waals gap and the larger cleavage energy[7,25,28,31], 


which may induce a 3D magnetic phase transition. Forthe purpose of confirmation, 
we performed a detailed characterization of the critical phenomena using the initial 


isothermal M(H) curves around Tc for the CrGeTes3, which are shown in Fig.2(a). 
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Figure 2. (Color online): (a) Initial magnetization of CrGeTe; around Tg; (b) Arrott plots of M? versus 

H/M (the M(H) curves are measured at temperature intervals of 1 K and 0.5 K approaching Tc);(c) 

Normalized slopes as a function of temperature; (d) Modified Arrott plot (M'’versus (H/M)"”) of 


isotherms withf =0.24 and y =1 for CrGeTe3. The red dashed line is the linear fit of isotherm at 67.9 K; 


72 7e 


(e) Temperature dependence of Msand 7, t. The Tc and critical exponents are obtained from the fitting 


of Eq. (S1) and (S2); (f) The Kouvel-Fisher plot. The Tc and critical exponents are obtained from the 
linear fit. 


In the mean-field theory, the critical exponents and critical temperature can be 
determined from the Arrott plot withf of 0.5 and y of 1.0[42,43]. According to this 
method, the M? versus H/M (shown in Fig.2(b)) should be a series of parallel straight 
lines in the higher field range around Tc and the line at T = Tc shouldpass through the 
origin. Note that the lower field data mainly represent the arrangement of magnetic 


domains, which should be excluded in the fitting process[44]. However, all the curves 
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in Fig.2 (b) show nonlinear behaviors having downward curvature even in high field, 
which indicates non-mean-field-like behavior. Moreover, the positive slope reveals a 
second-order phase transition according to the criterion proposed by Banerjee[45]. As 
such, a modified Arrott plot should be employed to obtain the critical exponents. 

To determine an accurate model, we obtain a modified Arrott plot following Eq. 
(S9) for single-crystallineCrGeTe; at different temperatures. Three groups of possible 
exponents belonging to the 3D Heisenberg model (6 = 0.365, y = 1.386), 3D-Ising 
model (£ = 0.325, y = 1.24) and tricritical mean-field model (4 = 0.25, y = 1.0) exhibit 
nearlystraight lines in the high field region[46,47]. We calculate their normalized 
slopes (NS)defined as NS = S(T) / S(T-""° =67.3K) .By comparing NS with the ideal 
value of unity, one can identify the most suitable model[46,47]. Figure 2(c)shows the 
plots of NS versus Temploying the three different models, revealing that the tricritical 
mean-field model is the most appropriate to describe the critical behavior of 
CrGeTe3.By proper selections of / and y, one can clearly show the isotherms are a set 
of parallel straight lines in high fields as displayed in the Fig.2(d). The linear 
extrapolation from the high field region gives the spontaneous magnetization M;(T,0) 
and the inverse of initial susceptibility 7,'(T,0) (seeFig.2 (e)) corresponding the 
interceptsonthe M"? and (H/M)"axes, respectively. By fittingthe data ofMs(T,0) and 
XY (T,0)to Eq. (S1) and (S2), one obtains two new values of £ = 0.242 + 0.006 with 
Tc = 67.95 + 0.01 and y = 0.985 + 0.009 with Tc = 67.90 + 0.09. These results are 
again very close to the critical exponents of tricritical mean-field model. In addition, 


these critical exponents and Tc can be obtained more accurately from the 


Kouvel-Fisher (KF) method[48]. Hence, onecan find that the temperature dependence 
of M,(dM,/dT) ‘and y,'(dy,'/dT) ‘should be straight lines with slopes 1/8 and 
1/y, respectively. As seen in Fig.2(f), the linear fit yieldsthefof 0.240 + 0.006 with 
Tcof 67.91 + 0.07 and yof 1.000 + 0.005 with Tcof 67.88 + 0.05, respectively. 
Remarkably, the obtained values of the critical exponents and Tc using the KF method 
are in excellent agreement with those using the modified Arrott plot based on 
thetricritical mean-field model. This suggests that the estimated values are 
self-consistent and unambiguous. 

To further validate the above critical exponents f andy, we study the relation 
among these exponents. According to Eq. (S3), 6 can be directly estimated from the 
critical isotherm at Tc. Figure 3 (a) shows the isothermal magnetization M(H) at Tc = 
67.9 K. The inset of the same plot has been demonstrated on a log-log scale. The 
solid straight line with a slope 1/6 is the fitted result using Eq. (S3). From the linear fit 
we obtained the third critical exponent dof 5.032 + 0.005.Moreover, the exponent ô 
can be calculated by the Widom scaling relation[49,50] 

d=1+7/B. (3) 
Based on the f# and y values calculated in Fig.2(e) and (f), Eq. (8)yields do0f5.070 + 
0.006 and 5.167 + 0.006, respectively. We emphasize that these values are very close 
to the results from the experimental critical isothermal analysis. Therefore, the critical 
exponents obtained in this study basically obey the Widom scaling relation, showing 
that the obtained /,yand dare reliable. 

Finally, these critical exponents should follow the scaling equation (Eq. (S6)) in 
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the critical region. The scaling equation indicates that m versus h forms two universal 
curves for T>Tc and T<T¢, respectively. Based on Eq. (S7), the isothermal 
magnetization around the critical temperatures for CrGeTe3 has been plotted in Fig.3 
(b). All experimental data in the higher field region collapse onto two universal curves, 
in agreement with the scaling theory. The inset of Fig.3 (b) shows the corresponding 
log-log plot. Similarly, all the points also collapse into two different curves in the 
higher field region. This result shows again that the obtained results of the critical 


exponents and Tc are valid. 


‘[b) 


Bk 


22 


0 
4 ô= 5.032(5) 


6 12 = 18 24 
H (kOe h = Hle{""” (10° Oe) 


Figure 3.(Color online): (a) Isothermal M(H) at Tc. The inset shows the alternative plot on a log-log 
scale and the straight line is the linear fit following Eq. (7); (b) Renormalized magnetization m versus 


renormalized field h at several typical temperatures around the Tc. The inset shows the alternative plot 
on a log-log scale; the effective exponents (c) Px below Tc and (d) 7 ,¢¢ above Tc as a function of the 


reduced temperature €. 


To further examine the convergence of the critical exponents, the effective 


exponents f, and 7,,, can be are obtained by Eq. (S8) and (S9) for CrGeTe3. As shown 
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in Fig.3 (c) and (d), both #,,, and y, show a non-monotonic variation with € (see Eq. 
(S4)). The lowest € (€min) are 5.89X10° and 1.47x10° for Sẹ and 7.4, respectively. 
We obtain the effective exponents Jẹ of 0.242 and y, of1.069, indicating that both 
Bas and yg are converged when the temperature approaches Tc. 

The experimentalcritical exponents of CrGeTe3, as well as the theoretical values 
of CrSiTes, MnSi and some other manganitesbased on various models, are 
summarized in Table 1. It is seen that the critical exponents for MnSi and doped 
manganites are consistent with those of tricritical mean-field theory[46,47,51,52]. 
These compounds have the same characteristics, i.e., a tricritical pointseparatingthe 
first-order from the second-order ferromagnetic phase transitions. This phenomenon 
shows that the element substitution[47,52], hole or electric doping[51], and external 
magnetic field[46] can induce the tricritical behavior. However, CrXTe3 presents a 
second-order ferromagnetic phase transitions[7,25-28] and our results indicate that the 
critical behavior of CrGeTe3 is close to the theoretical value of tricritical mean-field 
model. Comparing with CrSiTe3, showing characteristic of 2D-Ising model[7,34,35], 
the smaller van der Waals gap and the larger planar nearest-neighbor Cr-Cr distance of 
CrGeTezenhances the Curie temperature from 32K for the CrSiTe3 to 61K for the 
CrGeTe3[7,25-28]. In addition, the neutron scattering and isothermal magnetization 
experiments yield a critical exponent £ of around 0.151 or 0.17 for CrSiTe3[7,34,35], 


which is close to the value expected for a 2D transition ( 3}""* = 0.125) and well below 


the values expected for a3D transition (33° = 0.326).Our results yield a critical 


exponent p of 0.24 for CrGeTes that is close to the critical exponent of the tricritical 
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mean-field model. Hence, the increase of the X atom radius, facilitating super 


exchange coupling between the Cr atoms via the Te atom andleading to the smaller 


van der Waals gap in CrXTe3 system [7,25-28], could induce a tricriticalmagnetic 


phase transition in the CrGeTe3 single crystal. 


Table 1.Critical exponents of CrGeTe; with various theoretical models, CrSiTe3 and other related 


materials with tricritical mean-field model (SC = single crystal; PC = polycrystalline; cal = Calculated 


from Eq. (3)). 


Composition Reference | Tc (K) 
Technique Bp y ô 
CrGeTe;°~ This work | 67.9 Modified | 0.242 + | 0.985 + | 5.070 + 
Arrott plot | 0.006 | 0.003 | 0.006%" 
Kouvel-Fish | 0.240 + | 1.000 + | 5.167 + 
ermethod | 0.006 | 0.005 | 0.006%" 
Criticalisoth - - | 5.032 + 
erm 0.005 
Tricritical [45] - Theory 0.25 1 5 
mean-field 
Mean-field [43,44] - Theory 0.5 1 3 
3D-Heisenberg [43,44] - Theory 0.365 1.386 4.8 
theory 
3D-Ising [43,44] z Theory 0.325 1.24 4.82 
CrSiTe;°° [34] 31 Modified | 0.170 + | 1.532 + | 10.012 
Arrott plot 0.008 0.001 Ł 
0.047%" 
Mnsi*< [46] 30.5 Modified | 0.242 + | 0.915 + | 4.734 + 
Arrott plot 0.006 0.003 0.006 
Lao.1Ndo.¢Sto.3MnO3 [47] 249.3 | Modified | 0.257+)| 112+ | 5.17+ 
vs Arrott plot | 0.005 | 0.03 0.02 
Lao. 9Teo,;MnO;"° [51] 239.5 | Modified | 0.201+) 1.27+ | 7.14+ 
Arrott plot 0.003 0.04 0.04 
Lag ¢Cap 4MnO;"° [52] 265.5 | Modified | 0.25+ | 1.034 | 5.04 
Arrott plot 0.03 0.05 0.8 
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Although single-crystallineCrSnTe; has not yet been synthesized, we speculate 
that the magnetism of CrSnTe3 should be closer to the 3D-Ising model. To support this 
assumption, we perform DFT calculations with the same calculation parameters that 
were used in Ref. [53]. Figure 4 (a) shows the calculated formation energyF;, which is 
defined as the energy cost of extracting a sheet of single-layer CrXTe3 from their bulk 
counterparts.As can be seen, E; increases as the species vary from Si to Ge. This is 
consistent with the larger theoretical cleavage energy of single-layer CrGeTe3than that 
of CrSiTe3, which indicates that the layers are coupled more strongly in CrGeTe3[31]. 
The formation energy of CrSnTe3; is even higher than the other two 
compounds,revealing that it presentsthe strongest interlayer coupling, which leads to 
its 3D characteristics.Figure 4 (b), (c), and (d) illustrate the charge density of the three 
compounds. Consistent with the trend of the E;results, the electron density around the 
Sn-Sn pairis the least amongthe three materials. Namely, more electrons in 


CrSnTesparticipate the interlayer coupling. 
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Figure 4. (Color online): (a) Formation energy of single-layer CrXTe;; The formation energy of 
single-layer CrSiTe; is adopted from Ref.[53](b), (c), and (d) Charge density of bulk CrXTe3with an 


isosurface value of 0.05 e/rgon 


In conclusion, we have performed a comprehensive experimental study on the 
critical properties of single-crystalline CrGeTe3; using isothermal magnetization 
around the Curie temperatureT. Based on various experimental techniques including 
the modified Arrott plot, Kouvel-Fisher method and critical isotherm analysis, we 
obtained the critical exponents £, y, and 6 of 0.240 + 0.006, 1.000 + 0.005, and 5.070 
+ 0.006, respectively, at the Curie temperature of 67.9 K.These numerical results are 
similar to the theoretical values in the tricritical mean-field model, which is therefore 
capable of describing the critical magnetic behavior of 2D CrGeTe3.DFT calculations 
show that the formation energy of CrGeTe3 lies between those of CrSiTe3 and 
CrSnTe3, which isin line with a crossover of the magnetic phase transition from 2D to 
3D. Overall, our findings provide a fundamental understanding of the anomalous 
PM-FM transition in a novel 2D ferromagnetic semiconductor. 


Acknowledgements 


This work was supported by the National Key Research and Development Program 
under contracts 2016YFA0401803 and 2016YFA0300404, the Joint Funds of the 


15 


National Natural Science Foundation of China and the Chinese Academy of Sciences' 
Large-Scale Scientific Facility under contract U1432139and U1532153,the National 
Natural Science Foundation of China under contract11674326, Key Research 
Program of Frontier Sciences, CAS (QYZDB-SSW-SLH015), and the Nature Science 
Foundation of Anhui Province under contract 1508085ME103.J. Z. is supported by 
the Nature Science Foundation of China (Grant No. 51602079) and the Fundamental 
Research Funds for the Central Universities of China (Grant No. 372 
AUGA5710013115). This research also used computational resources of the National 
Supercomputing Center of China in Shenzhen (Shenzhen Cloud Computing Center). 


References 

[1] K. S. Novoselov, A. K. Geim, S. V. Morozov, D. Jiang, Y. Zhang, S. V. Dubonos, I. 
V. Grigorieva, and A. A. Firsov, Science 306, 666 (2004). 

[2] K. S. Novoselov, A. K. Geim, S. V. Morozov, D. Jiang, M. I. Katsnelson, I. V. 
Grigorieva, S. V. Dubonos, and A. A. Firsov, Nature 438, 197 (2005). 

[3] Y. Zhang, Y.-W. Tan, H. L. Stormer, and P. Kim, Nature 438, 201 (2005). 

[4] A. K. Geim and K. S. Novoselov, Nat. Mater.6, 183 (2007). 

[5] A. H. Castro Neto, F. Guinea, N. M. R. Peres, K. S. Novoselov, and A. K. Geim, 
Rev. Mod. Phys. 81, 109 (2009). 

[6] N. Zibouche, A. Kuc, J. Musfeldt, and T. Heine, Ann. Phys. (Berlin) 526, 395 
(2014). 

[7] T. J. Williams, A. A. Aczel, M. D. Lumsden, S. E. Nagler, M. B. Stone, J. Q. Yan, 
and D. Mandrus, Phys. Rev. B 92, 144404 (2015). 

[8] L Žutić, J. Fabian, and S. Das Sarma, Rev. Mod. Phys. 76, 323 (2004). 

[9] A. H. MacDonald, P. Schiffer, and N. Samarth, Nat. Mater. 4, 195 (2005). 

[10] T. Dietl, Nat. Mater. 9, 965 (2010). 

[11] Z. Deng, C. Q. Jin, Q. Q. Liu, X. C. Wang, J. L. Zhu, S. M. Feng, L. C. Chen, R. 
C. Yu, C. Arguello, T. Goko, F. Ning, J. Zhang, Y. Wang, A. A. Aczel, T. Munsie, T. J. 
Williams, G. M. Luke, T. Kakeshita, S. Uchida, W. Higemoto, T. U. Ito, B. Gu, S. 
Maekawa, G. D. Morris, and Y. J. Uemura, Nat. Commun. 2, 422 (2011). 

[12] Z. Y. Zhu, Y. C. Cheng, and U. Schwingenschlögl, Phys. Rev. B 84, 153402 
(2011). 

[13] D. Xiao, G.-B. Liu, W. Feng, X. Xu, and W. Yao, Phys. Rev. Lett. 108, 196802 
(2012). 

[14] W. Feng, Y. Yao, W. Zhu, J. Zhou, W. Yao, and D. Xiao, Phys. Rev. B 86, 165108 
(2012). 

[15] W.-Y. Shan, H.-Z. Lu, and D. Xiao, Phys. Rev. B 88, 125301 (2013). 


16 


[16] X. Xu, W. Yao, D. Xiao, and T. F. Heinz, Nat. Phys. 10, 343 (2014). 

[17] W. Han, R. K. Kawakami, M. Gmitra, and J. Fabian, Nat. Nano. 9, 794 (2014). 
[18] H. Ohno, D. Chiba, F. Matsukura, T. Omiya, E. Abe, T. Dietl, Y. Ohno, and K. 
Ohtani, Nature 408, 944 (2000). 

[19] C.-Z. Chang, J. Zhang, X. Feng, J. Shen, Z. Zhang, M. Guo, K. Li, Y. Ou, P. Wei, 
L.-L. Wang, Z.-Q. Ji, Y. Feng, S. Ji, X. Chen, J. Jia, X. Dai, Z. Fang, S.-C. Zhang, K. 
He, Y. Wang, L. Lu, X.-C. Ma, and Q.-K. Xue, Science 340, 167 (2013). 

[20] X. Kou, S.-T. Guo, Y. Fan, L. Pan, M. Lang, Y. Jiang, Q. Shao, T. Nie, K. Murata, 
J. Tang, Y. Wang, L. He, T.-K. Lee, W.-L. Lee, and K. L. Wang, Phys. Rev. Lett. 113, 
137201 (2014). 

[21] Z. Wang, C. Tang, R. Sachs, Y. Barlas, and J. Shi, Phys. Rev. Lett. 114, 016603 
(2015). 

[22] H. Gonzalez-Herrero, J. M. Gomez-Rodriguez, P. Mallet, M. Moaied, J. J. 
Palacios, C. Salgado, M. M. Ugeda, J.-Y. Veuillen, F. Yndurain, and I. Brihuega, 
Science 352, 437 (2016). 

[23] C. Gong, L. Li, Z. Li, H. Ji, A. Stern, Y. Xia, T. Cao, W. Bao, C. Wang, and Y. 
Wang, Nature (2017). 

[24] N. Sivadas, M. W. Daniels, R. H. Swendsen, S. Okamoto, and D. Xiao, Phys. Rev. 
B 91, 235425 (2015). 

[25] V. Carteaux, D. Brunet, G. Ouvrard, and G. Andre, J. Phys.: Condens. Matter 7, 
69 (1995). 

[26] H. Ji, R. A. Stokes, L. D. Alegria, E. C. Blomberg, M. A. Tanatar, A. Reijnders, L. 
M. Schoop, T. Liang, R. Prozorov, K. S. Burch, N. P. Ong, J. R. Petta, and R. J. Cava, 
J. Appl. Phys. 114, 114907 (2013). 

[27] Y. Tian, M. J. Gray, H. Ji, R. Cava, and K. S. Burch, 2D Materials 3, 025035 
(2016). 

[28] H. L. Zhuang, Y. Xie, P. R. C. Kent, and P. Ganesh, Phys. Rev. B 92, 035407 
(2015). 

[29] B. Siberchicot, S. Jobic, V. Carteaux, P. Gressier, and G. Ouvrard, J. Phys. Chem. 
100, 5863 (1996). 

[30] S. Lebégue, T. Björkman, M. Klintenberg, R. M. Nieminen, and O. Eriksson, 
Phys. Rev. X 3, 031002 (2013). 

[31] X. Li and J. Yang, J. Mater. Chem. C 2, 7071 (2014). 

[32] J. Zhang, B. Zhao, Y. Yao, and Z. Yang, Phys. Rev. B 92, 165418 (2015). 

[33] J. Liu, S. Y. Park, K. F. Garrity, and D. Vanderbilt, Phys. Rev. Lett. 117, 257201 
(2016). 

[34] B. Liu, Y. Zou, L. Zhang, S. Zhou, Z. Wang, W. Wang, Z. Qu, and Y. Zhang, Sci. 
Rep. 6, 33873 (2016). 

[35] V. Carteaux, F. Moussa, and M. Spiesser, EPL Europhys. Lett. 29, 251 (1995). 
[36] G. Kresse and J. Furthmüller, Phys. Rev. B 54, 11169 (1996). 

[37] D. M. Ceperley and B. J. Alder, Phys. Rev. Lett. 45, 566 (1980). 

[38] J. P. Perdew and A. Zunger, Phys. Rev. B 23, 5048 (1981). 

[39] P. E. Blöchl, Phys. Rev. B 50, 17953 (1994). 

[40] G. Kresse and D. Joubert, Phys. Rev. B 59, 1758 (1999). 


17 


[41] C. Kittel, Introduction to Solid State Physics. (Wiley, 2004). 

[42] A. Arrott, Phys. Rev. 108, 1394 (1957). 

[43] S. N. Kaul, J. Magn. Magn. Mater. 53, 5 (1985). 

[44] N. Khan, P. Mandal, K. Mydeen, and D. Prabhakaran, Phys. Rev. B 85, 214419 
(2012). 

[45] B. K. Banerjee, Phys. Lett. 12, 16 (1964). 

[46] L. Zhang, D. Menzel, C. Jin, H. Du, M. Ge, C. Zhang, L. Pi, M. Tian, and Y. 
Zhang, Phys. Rev. B 91, 024403 (2015). 

[47] J. Fan, L. Ling, B. Hong, L. Zhang, L. Pi, and Y. Zhang, Phys. Rev. B 81, 144426 
(2010). 

[48] J. S. Kouvel and M. E. Fisher, Phys. Rev. 136, A1626 (1964). 

[49] B. Widom, J. Chem. Phys. 43, 3898 (1965). 

[50] B. Widom, J. Chem. Phys. 41, 1633 (1964). 

[51] J. Yang, Y. Lee, and Y. Li, Phys. Rev. B 76, 054442 (2007). 

[52] D. Kim, B. Revaz, B. L. Zink, F. Hellman, J. J. Rhyne, and J. F. Mitchell, Phys. 
Rev. Lett. 89, 227202 (2002). 

[53] M.-W. Lin, H. L. Zhuang, J. Yan, T. Z. Ward, A. A. Puretzky, C. M. Rouleau, Z. 
Gai, L. Liang, V. Meunier, B. G. Sumpter, P. Ganesh, P. R. C. Kent, D. B. Geohegan, 
D. G. Mandrus, and K. Xiao, J. Mater. Chem. C 4, 315 (2016). 


18 


